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r^ Abstract 

fj We report measurements of the phonon density of states as measured with inelastic x-ray scattering in SmFeAsOi-jF, powders. An 
I unexpected strong renormalization of phonon branches around 23 meV is observed as fluorine is substituted for oxygen. Phonon 
f^ dispersion measurements on SmFeAsOi_iFy single crystals allow us to identify the 21 meV Aig in-phase (Sm,As) and the 26 meV 
^ Big (Fe,0) modes to be responsible for this renormalization, and may reveal unusual electron-phonon coupling through the spin 
^ channel in iron-based superconductors. 
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' The recent discovery of superconductivity in FeAs based 
compounds has sparked a new gold rush amongst the strongly 
correlated electron community yj]. Among the various fami- 
lies of FeAs based superconductors, the so-called '1111', with 
composition REFeAsO (RE being a rare earth), has the high- 
est values of the superconducting transition temperatures Tc'. 
depending on the lanthanide, Tc can be as high as 55 K upon 
fluorine doping \^. The high value of Tc, early DFT calcu- 
lations yD, and other physical properties, appeared to support 
jthe notion of an exotic mechanism leading superconductivity 
■in these materials. There is however increasing speculation that 
an exotic electron-phonon coupling, possibly enhanced through 
the spin or orbital channels, is at play in the iron-based high- 
temperature superconductors. In particular, the fact that the Fe- 
As configuration is intimately related to Tc ISO and to the mag- 
netic state of the iron sublattice la |7[] show the lattice is active 
in this problem. This was further confirmed recently through 
the observation of strong renormalization of zone center Ra- 
man 18i|9|] or infrared jlOll active phonons at Tai in BaFe2As2. 
Here, we present inelastic x-ray scattering (IXS) measurements 
on poly- and single-crystalline samples of SmFeAsOi-jF, that 
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allowed us to investigate the doping dependence of the gener- 
alized phonon density of states (GPDOS) and of the phonon 
dispersions. We report an unexpectedly strong doping-induced 
renormalization of two phonon branches - the 2 1 me V ( As,Sm) 
and 26 meV (Fe,0) modes - providing significant insight into 
the evolution with doping of the energy and momentum depen- 
dence of the coupling between the lattice and the spin degrees 
of freedom, and suggesting that the spin-phonon coupling may 
be of key importance in these compounds. 

2. Results 

SmFeAsCGo.gFo.i) samples (Sm-1111) was prepared under 
high pressure and high temperature using a "belt" type high 
pressure cell in a similar way to that reported previously for 
LaFeAsOi_iFv lllll . Sm, Fe, Fe203, As and SmFs (in the case 
of the fluorine doped sample) powders were mixed together and 
pressed in the form of cylindrical pellets. They were then intro- 
duced into a crucible (machined from rods of hexagonal boron 
nitride) surrounded by a cylindrical graphite resistive heater and 
the whole assembly was placed in a pyrophyllite gasket. The 
samples were treated at 6GPa, 1000-1 100 °C for 4 hours, then 
quenched to room temperature and depressurized. The XRD 
patterns show that the major phase is the one expected, Sm- 
1111 with some small impurities of FeAs, SmAs and Sm203. 
SmFeAsO synthesis is described in ref. 11211 . 
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Figure 1 : a) comparison between the measured generalized phonon density of 
states of SmFeAsO (black curve) and SmFeAsOo.gFo.i (red curve), measured at 
room temperature. The area of both spectra have been normalized to unity. The 
blue line corresponds to the difference between the two spectra, b) partial PDOS 
obtained from DFT calculation Q] and weighted according to the procedure 
described in ref. 11611. 



The IXS experiment on polycrytalline SmFeAsO and 
SmPeAsOo.gFo.i was carried out on beamline ID28 at the ESRF 
using the silicon (9 9 9) configuration at 17.794 keV with an in- 
strumentation energy resolution of 3.0 meV. Within the limit of 
single phonon scattering (Stokes process), the intensity of the 
inelastically scattered signal with energy (momentum) transfer 
(^ {Q - {Qa QbQc)) is directly proportional to the dynamical 
structure factor lll3ll : 



SiilM = 2 X 



/;,(Q)e'Qr«-«'"(Q.,r„(Q,;)) 



(1 + n(u)Qj, T))6{a) - ujqj) 



'^Q. 



(1) 



where /„(Q) is the atomic form factor of the n''' atom (at posi- 
tion r„) in the unit cell, M„ its mass and W„ the Debye-Waller 
factor The summation index j refers to the /'' phonon branch 
and a"„(Q, j) and ojq to the associated eigenvector and eigen- 
values, respectively. Finally, «(w, T) - {exp{cL)/kBT) - 1)"' 
stands for the Bose factor The scattered photons are analyzed 
by a set of 8 analyzers that were centered on two different set 
of angles to cover a Q range from 53 nm"' to 73 nm"' ({Q) = 
63 nmr'^), and in this case the dynamical structure factor must 
be averaged over th e sp here of radius Q - ||Q||. As previ- 
ously demonstrated pi 41], averaging several measurements over 



a large range of Q {Q„,i„ <Q < Q,„ax) gives access to: 

S(ciJ) oc xGPDOS(cj), (2) 



with the GPDOS being defined as a weighted sum of the par- 
tial phonon densities of states, G„(w), for each atomic species 



GPDOS (lo) = Yj 



M„ 



(3) 



In the upper panel of Figure [T] we show the quantity 
Imeiasiic{<^) 'x[a)/(l+ n(u), T))] measured at room temperature on 
SmFeAsO and nominally doped SmFeAsOo.gFo.i. The inelastic 
signal Imeiasiici(^) has been obtained after summation of the indi- 
vidual spectra weighted by the appropriate analyzer efficiency, 
subtraction of the elastic contribution, and recursive correc- 
tion from multiphonon contributions, following refs. il4 Il5n . 
As previously reported for the case of NdFeAsOi_;tFx UoH 
and LaFeAsOi-.iFv U7\, we observe essentially three peaks of 
about 5-6 meV full-width-half-maximum located around 13, 
24 and 32 meV. In the lower panel of Figure 11] we report 
the theoretical results obtained by Boeri et al. [3'] on a sim- 
ilar compound, LaOFeAs. To allow a direct comparison be- 
tween the experimental data and the theoretical calculation, we 
have weighted each calculated partial PDOS G„{co) by its cor- 
responding e''^^" fniiQYf' M„^ factor (these factors are given in 
ref. il6ll ). The result of this procedure has been plotted in the 
lower panel of Fig. [1] The overall agreement with the experi- 
mental data for the parent compound is good for the two first 
peaks, but as previously observed in the cases of NdPeAsO and 
SmFeAsO, the higher energy peak is found to be harder, by 
about 4 meV, calculated than is observed in the experiment. 

The GPDOS measured on SmFeAsOo.gFo.i and plotted in the 
upper panel of Figure [1] reveal a spectra with a similar overall 
shape to the one measured on SmFeAsO, except for the second 
peak. There is indeed, a striking change below 24 meV, with a 
clear increase of the spectral weight around 21 meV similar to 
our previous findings in NdFeAsOi_vFv lll6ll . 

According to first principle calculations, this middle band 
originates from three optical phonons branches extending from 
20 to 28 meV. These modes are Raman active at the zone cen- 
ter and are associated with the in-phase A\g (As,Sm) vibration 
(~ 21 meV at Q = in SmFeAsO from Raman measure- 
ments lll8l 11911 ). the out-of-phase Aij (As,Sm) mode (Ai^(2) 



- 23 meV) and the Bi^ (Fe,0) mode (26 meV), respectively. 

To go further and determine which of these modes are 
responsible for the observed renormalization, we have car- 
ried out IXS experiments on SmFeAsO and superconducting 
SmFeAsOo.6Fo.35 single crystals to measure their phonon dis- 
persion, with particular emphasis on these three modes. High 
quality single crystals of SmFeAsO and SmFeAsO0.60F0.35 
of about 100x100x20 yum-' were grown as described else- 
where II20II21I1 . Our SQUID measurements on SmFeAsO show 
magnetic ordering at Tt^ =130 K, and the superconducting 
transition in the nominally doped SmFeAsO0.60F0.35 takes place 
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Figure 2: Experimental IXS spectra of SmFeAsOi_,Fy for q =(0.15 0.15) 
(left) and q =(0.35 0.17) (right) measured in the Brillouin zones centered 
around F = (1 12). The red lines are the result of the total fit, and the thin gray 
lines the individual phonons. The in-phase (red) and out-of-phase (green) Aig 
(As,Sm) modes as well as the Big (Fe,0) mode (blue) have been highlighted. 
The black lines represent the theoretical spectra calculated for LaFeAsO at 
these Q vectors from the DFT 12411 . 



at Tc - 53 K. We used an incident photon energy of 21.747 
keV (silicon (11 11 11) configuration) providing an energy res- 
olution of 2.0 meV on ESRF beamline ID28. The x-ray beam 
was focused down to 50x40 //w". Rocking curves of 0.5° were 
measured on the (0 4) reflection, indicating a good c-axis mo- 
saicity. 

From eq. [1] one can see that probing c-axis polarized 
phonons, such as of those we are interested in, requires one 
to carry out measurements in a zone with the largest possible 
Qc component along the c* reciprocal lattice vector (note that 
in what follows, Q is defined with respect to the tetragonal unit 
cell, and q - (qa qb qc) refers to the reduced momentum trans- 
fer relative to the zone center). To enhance the contrast of the 
optical modes with respect to the strong acoustic one, we also 
decided to work close to a weak Bragg spot. The best com- 
promise was found to be the zone around F = (1 12). This 
corresponds to a challenging quasi-grazing incidence reflection 
geometry (because of the size of the crystals), with a scattering 
angle 20 ~ 49°. This angle is close to the highest achievable 
for 1D28, and we had to perform the transverse scans along the 
100 direction. 

In Fig. |2] we present typical IXS spectra obtained for Sm- 
FeAsO and SmFeAsOo.esFojs, together with the result from 
DFT calculation performed for LaFeAsO from ref . ll24ll . Exper- 
imental data were analyzed by fitting to a series of Lorentzians, 
as shown in Fig. |2] With the exception of the elastic line, fit- 
ting of the data to resolution limited peaks gave rather poor re- 
sults, and slightly broader lineshapes (2.2 to 2.5 meV) were 
used. Similar broadening has also been reported in the case of 
CaFe2As2 112211 . The assignment of the phonon branches is chal- 
lenging due to the large discrepancies between the IXS inten- 
sities predicted by the DFT calculations and the experimental 
spectra, as illustrated in Fig. |2] We by-passed this difficulty, 
at least for the three c-axis polarized branches around 23 meV 
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Figure 3: Momentum dependence of the doping induced renormalization of the 
26 meV, Bij (upper panel) and 21 meV, in-phase Aig (lower panel) modes in 
the {qa ijr) plane. The empty circles correspond to the analyzer positions, the 
dashed ellipse indicates the q-resolution, and the white areas where there is in- 
sufficient data. The insets show schematics of the respective q=0 eigenvectors, 
with the atoms color coded as indicated (from ref. I23II V 



we are interested in, by taking advantage of measurements per- 
formed in an adjacent Brillouin zone (centered around F = 1 
1 1), as described in ref. |23]. 

Using this method and exploiting the 8 analyzers available on 
ID28, we have been able to map out the dispersion of the three 
c-axis polarized modes in the {q^, qc) plane for the parent and 
doped compounds, at room temperature i23n . The renormal- 
ization of the modes was found to be rather weak in the FX 
direction, but displayed a surprisingly strong dependence with 
qc- As shown in Fig.|2l we observed a softening of the in-phase 
Aig mode and a hardening of the Big mode. The maximum am- 
plitude for these effects are obtained close to q=(0.3 0.3) for 
the in-phase Ai^ mode (1.2 meV softening) and close to q=(0.5 
0.5) for the Big mode (1.7 meV hardening). No clear doping 
dependence for the out-of-phase Aig (As,Sm) mode was found. 

3. Discussion 

The disagreement between the IXS intensities predicted by 
standard DFT calculation and experimental data has also been 
reported in the case of the 122 compounds 1221 |25L 12611 . In 
these cases a much better agreement between theory and exper- 
iment in parent compounds above their magnetic transition has 
been found by including spin-polarization in the calculations. 
Ref. I26n also mentions the fact that it improves the agreement 
between the experimental and calculated GPDOS in the case of 
the 1111 compounds. This indicates that magnetism is the miss- 



ing ingredient in the ab-initio calculations discussed above. It 
was also found that the magnetic ground state was intimately re- 
lated to the Fe-As distance along the c-axis f6', 7], which would 
naturally lead to a coupling between the atomic motion mod- 
ulating this distance - precisely the in-phase out-of-plane Ai,, 
(As,Sm) and the (Fe,0) Bi^ modes - and the amplitude of the 
Fe moments 112711 . 

As far as the doping dependence of the phonon spectra is 
concerned, a calculation based on the virtual crystal approxima- 
tion suggested that doping should induce only minor changes 
in the phonon spectra yD. Another study taking into account 
the F-doping explicitly in a 2x2x1 supercell within an LDA ap- 
proach |24|l has related the doping-induced changes seen in the 
PDOS lll6ll to structural relaxation (with a relative change of 
the unit cell volume AV/V of about 4 %). However, this ap- 
proach leads to a hardening of the in-phase (As,Sm) Ai^ mode 
that contrasts with the observed softening, and overestimates 
the structural changes as experimentally, the doping-induced 
unit cell volume is rather limited {AV/V ~ 0.5%). The ob- 
served doping-induced phonon renormalization and their mo- 
mentum dependencies reported here may rather reflect the evo- 
lution upon fluorine doping of momentum and frequency of ei- 
ther the electron-phonon coupling constant g(q, aS) or the elec- 
tronic/magnetic susceptibility xi^l^ (^) of the system. To fully 
understand the role played by the coupling of phonon to mag- 
netism in iron-pnictides, further theoretical and experimental 
investigations are needed. We also note that in a recent calcula- 
tion for the case of 122 compounds, it was shown that, in addi- 
tion to a renormalization of the phonon frequencies, magnetism 
also tends to enhance the electron-phonon coupling constant, 
though this does not appear to be sufficient to explain the high 
superconducting Tc of these materials 1281]. 



4. Conclusion 

In summary, we have found a clear doping dependence 
for two c-axis polarized phonon branches in the 1111 iron- 
pnictide family that may reveal unusual electron-phonon cou- 
pling through the spin channel in iron-based superconductors. 
Further theoretical and experimental investigations are needed 
to fully understand the role played by such coupling in these 
materials properties. 
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